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^ ' Abstract 

-)-* 

The ac penetration depth, Xac{T, H, /), has been measured in superconduct- 

, . ing UPt3 single crystals for H\\c and ii"_Lc in the range / =0.01-1 MHz and 

j-H , T =0.1-0.4 K. The contributions from bulk and surface pinning have been 

O I separated to yield the flux-flow resistivity Pi{H). With H ± c, pf displays 

magneto-resistance at low-field which agrees with previous dc measurements 

and the characteristic scaling law of clean crystals with anisotropic gap. When 

^ . -'^ II c, the low-field pf is three times larger. We interpret this property as ev- 

5^ i idence for flux lines with unconventional core structure. 
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The reduced symmetry of the superconducting energy gap in unconventional supercon- 
ductors is generally inferred from the presence of low-energy excitations, which exist due to 
nodes in the gap function A{kp). In the most important heavy-fermion compound UPta, ev- 
idence for gap nodes has been deduced from measurements of the temperature dependences 
of the London penetration depth A(T) |1|, ultrasonic attenuation 0, thermal conductivity 
0], and most recently flux- flow resistivity Pi{T) 0. UPts also has a unique phase diagram, 
with two low- field phases A and B and a high-field C-phase [^Q (see insets in Fig. ^. These 
properties are compatible with a Ginzburg-Landau expansion in terms of a two-component 
order parameter |0-0. Such states allow for new and unconventional vortex core structures 
as well as other topological defects, such as domain walls or textures |7|-p!0|j. A multitude of 
different vortex structures and other topological defects have been discovered in ^He super- 
fluids [^, but so far direct experimental evidence for their existence in superconductors is 
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sparse 

From neutron diffraction measurements it is known that the B and C phases in UPta 
support distorted lattices of singly quantized vortex lines, with a clear break in the char- 
acteristic length scales at the B —>■ C transition |T^. At low magnetic field parallel to the 
c-axis, the two-component order-parameter model predicts that three different core struc- 
tures of singly quantized vortices are possible in the highly symmetric i?-phase [0. These 
have all singular structure with a radius comparable to the coherence length ^. One of them 
is axisymmetric while the other two, called triangular and crescent cores, are nonaxisym- 
metric and in contrast to the former, here the order parameter amplitude does not go to 
zero anywhere within the core. The relative stabilities of these different structures depend 
on the Ginzburg-Landau parameters, which are still largely undetermined. 

The investigation of the vortex-core structure in the bulk superconductor requires mea- 
surement of dynamical properties, like flux-flow dissipation or the Hall angle. We report on 
measurements of flux-flow resistivity pf{H,T) in the low-temperature limit (T ^ 0.3 K) for 
H A- c and H \\ c. For H J^ c our results are consistent with those of Ref. [^] at T ^ 0.3 
K. When H \\ c, the resistivity proves to be much larger than expected. It does not appear 
possible to explain this anisotropy unless the vortex structure is different in the two crystal 
directions. 

Vortex-dynamic measurements are complicated by pinning due to sample defects, which 
strongly alter the dc and low-frequency ac responses [|T^. At low temperatures in the B- 



phase, the use of large dc currents, to depin vortex lines and to bias the system into the 
free-flow state, is prohibited due to excessive Joule heating. In the present ac measurements 
this difficulty is overcome by recording the linear response over two decades in frequency 
(/ = 0.01-1 MHz) while the dissipation is maintained below 0.1 nW/mm^. As shown in 
Ref. |l^ , the linear regime at low rf excitation level (^ 1 /iT) and small amplitudes of 



vortex line oscillation (^ 1 A) lends itself to quantitative analysis, in spite of strong surface 
and/or bulk pinning. From the measurement of the effective ac penetration depth Aac(/) 
we extract the damping parameter, the flux-flow resistivity pf or the free-flow skin depth 



S{ = v/pf/(7r/io/). The frequencies are still sufficiently low so that the anomalous skin-depth 
{5{ > lOyum^ /m ~ 5000 A) and relaxation-time effects (r < 0.2 ns [H) can be neglected. 

Our measurements have been performed on two large single crystals (labelled B22 and 
B3b). They were prepared in Grenoble and were spark cut from the same ingot. Their 
quality is very similar to that of the sample used in the earlier dc measurements of Ref. [ffl. 



After cutting, the samples were annealed, but no surface polishing was applied. Their 
final dimensions are: L{x) x W{y) x d{z) = 5.5(a*) x 2.9(a) x 1.16(c) mm^ (B3b) and 
5.5(c) X 3.04(a*) x 0.63(d) mm^ (B22). They have low residual resistivity, pn = 0.52 + 
1.44T^ + 0.02 (/io-ff)^ iiQcm {^qH in Tesla), for currents J _L c, which is the case in our 
measurements. In these and similar samples 0,§], the resistivity for J || c is smaller by a 
factor 0.33 - 0.37 (for T ^ Tc). The measured pn corresponds to a mean-free path l^ > 10,^, 
proving that the samples represent the moderately clean limit |TB|. 



The dc field ^qH < 3 T and the vortex lines are along the z direction. The excitation 
field, /ie~'^'^^*, is applied along the x direction so that the vortices oscillate in the xz plane, 
and currents and electric fields E are induced in the y direction. Due to the Lorentz-force 
anisotropy, Ey is concentrated on the two main faces z = ±d/2 so that the sample mimics 
the 1-dim geometry of an infinite slab. The measured signal is the flux $ac through a 15-turn 
pick-up coil wound directly around the sample in the x direction. The apparent complex 
penetration depth is defined as Aac = A' + iA" = [$ac(-f^)-$ac(0)]/(2/xo/iW^)- 

The sample with its pick-up coil and the slightly larger excitation solenoid are placed 
inside the mixing chamber of a '^He-'^He dilution refrigerator. The measurements are per- 
formed in the temperature range T = 0.1-0.4 K, by recording the spectrum Aac(/) at fixed 
field H and by moving from one field value to the next, while the temperature is kept con- 
stant with a feed-back loop. The temperature is monitored with a calibrated Ge resistance 
thermometer inside the mixing chamber. Its field dependence is adjusted by comparing the 
measured Hc2 values with those in Ref . . The amplitude and phase of the signal voltage 
are calibrated against the difference between the responses in the normal and the Meissner 
states. This normalization procedure yields a phase accuracy better than 1° and a resolution 
S\„r ^ 1 Aim [ITSl . 



Typical spectra of Xadf) are shown in Fig. [^. The complete data set consists of more 
than 200 different spectra. The real and imaginary parts of Aac can be accurately fit with 
the following formula []T3|,|T7| 
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Ls{H,T) and Xc{H,T) are two frequency-independent lengths, which describe surface and 
bulk pinning respectively. The high-frequency limit, Aac(/— '•oo) = (H-i)5f/2, corresponds to 
the ideal fiux-fiow response, while the low-frequency limit, the quasistatic response, Aac(/— »• 
0) = A'(0) = (A^^ + -Lg^)~^, is purely inductive and cannot discriminate between surface 
and bulk pinning. The relative weight of surface and bulk pinning becomes apparent in the 
crossover regime. 

The excitation field penetrates as the sum of two modes. The first mode, localized near 
the surface, is associated with strong screening currents, whose amplitudes are governed by 
surface roughness. If the bulk sample is free of defects, the screened field penetrates further 
over the free-flow depth 6{. This situation (Ac— ^-cxd) has been systematically the case in the 
conventional superconductors, which so far have been measured JISl. In contrast, if there 



are bulk defects, such as those usually introduced in classical theories of pinning |1J] , the 
bulk mode is strongly attenuated, but penetrates at low frequencies over a Campbell length 
Ac (<C(5f at low frequencies). 

Curiously enough, our moderately clean UPts crystals display a large bulk pinning 



strength 1/Ac, as evident from Fig. |^. The data in Fig. |] show that in both crystal orien- 
tations at low fields pinning is strong and dominated by the surface process. This is not 
surprising, since after spark cutting the surfaces are visibly rough. With increasing field, 
surface pinning falls down and finally vanishes at a reversible field value close below the 
B ^C transition. This result is consistent with recent observations in point contact spec- 



troscopy ||T8H20|, which suggest a suppression of the order parameter at the surface at higher 
fields. Bulk pinning decreases more slowly as a function of [Hc2 — H) and is the dominant 
source of pinning in the C phase. No clear anomaly can be distinguished at the B —^ C 
transition jT^. Nor do we observe hysteresis as a function of the field-sweep direction. Low- 
field vortex-creep measurements |^ have suggested that new mechanisms could be present, 
such as intrinsic pinning by domain walls in the bulk . If this is the source for the bulk 
pinning, then the domain walls have to persist in the C phase. Additional work, including 
variations in surface treatment, should provide better insight in the pinning processes. 

The normalized fiux-fiow resistivity, Pf/pn, is shown in Fig. ^. As expected for the 
moderately clean superconductor, Pf{H) exceeds the "normal-core limit", pi = p^H/ Hc2, in 
both crystal directions. When the contributions from bulk and surface pinning have been 
removed, the pi{H) data become sample independent, as it should be. In the direction H ±c, 
where comparison is possible, pf{H) agrees quantitatively with the earlier dc measurements 
of Ref. 0]. It is usual to define the initial slope in Fig. ^ as 

At low-field the scaling factor ro(T) obtains, in the perpendicular direction, the value ro_L = 
1.6 ± 0.15 in the low-temperature limit. In fact, as seen in Fig. ^, below 0.3 K^ 0.6 Tc, tq 
is already temperature independent. This is in agreement with the results in Ref. |^, when 
they are extrapolated to low temperatures, and testifies for good general consistency between 
two different measuring methods and samples. It should be noted that the temperature 
independence below 0.3 K also embraces the bulk and surface pinning strengths in Fig. ^. 

In effect, Eq. (^ states that, at low field, fiux-fiow dissipation is additive and therefore 
proportional to the vortex density, n = poH/ipQ. Assuming that we can identify (po/poHc2 
with 27r6^, where b is an effective core radius, as is the case for conventional Abrikosov 
vortices, then Eq. (H) can be written as 



Pn{H) 



r*27mb' . (3) 



Here r* = r^ for Abrikosov vortices. More generally, when Hc2 is not an appropriate scal- 
ing parameter due to Pauli paramagnetic limitation of Hc2 and/or possible changes in the 
core structure, r*{T) can still be considered a dimensionless factor which depends on the 
mechanism of dissipation. 

In the dirty limit (/m *^ o^^ has Tq ^ 1 and Eq. (0) reads Pf/pn ~ H/Hc2 (normal-core 

model). In the clean limit dissipation is governed by the relaxation of the quasiparticle 

excitations localized in the vortex core. We shall make use of the theory by Kopnin et al. 

2^. They take 6 = ^ as the size of the confining vortex-core potential, allow for anisotropy 



in ^{kp) over a spherical Fermi surface, and predict that 



r*(T) ~ a ^^ a = 3^"^^" (4) 

A^ax(T) ' ((l-(A;-^)2)A2(fc)) ■ ^^ 

The () brackets denote a weighted average over the Fermi surface. For an isotropic gap 
(s-wave superconductor) the dimensionless factor a is unity. If A(fc) has nodal structure, 
or at least strong anisotropy, then a is expected to be larger than unity and its value 
depends on the orientation of the vortex lines with respect to the crystal axes. If we take 



Amax(0) = 1.9fcB7'c (and r* = ro), Eq. (|) yields a_L = 1.6x 1.9~3.0, which agrees within 
the combined experimental precisions with the value 3.2 obtained in Ref. @]. The large a 
is as expected because of gap anisotropy. 

The analysis of the measurements in terms of Eq. (^ applies for conventional supercon- 
ductors with To = r*, or equivalently, 6^ = ^a(,c = V^o/(2vr/io-ffc2±)- Here -^£211 is not a good 
scaling parameter because of the Pauli paramagnetic limitation at low temperature in this 
direction. To compare the results in the two orientations, we use Eq. (|^) by taking b\=^a^c 
and 6m = ^^. Furthermore, taking the high-temperature anisotropy C,a = 0.6^c, we find the 
low-temperature values of r* to be r]_ = 1.6 and rji =4.7/0.6 = 7.8. Eq. @ gives that these 
correspond to a± = 3.0 and a\\ = 7.8 x 1.9 = 14.8, which represents a large anisotropy of 
a||/a^~5. 

An anisotropy in a as large as this cannot be explained by the UPts structure, by taking 
into account the possible anisotropies from the D^h point group symmetry in the expression 
of a in Eq. (^. For instance, let us consider the gap structure in the 2-dim representations 
Eig {A r^ {kzkx:, kzky}) and E2u {^^ {kzi^l — ky),2kzkxky}), which are appropriate to the 
A and B phases ^. For Eig {Ar^kz{kx±iky), _B-phase) we find ay = 1.25 a^ = 4.4. For E2u, 
which is strongly supported by experiment PJ23[|, there is no anisotropy at all: a\\ =a_L = 3.9. 
Thus the small jump of about 10% in pf/pn, which was observed in Ref. |Q at the A^B 
transition for vortices in the aa* plane, might be ascribed directly to anisotropy in crystal 
structure and symmetry breaking, but not the large value of roy/roj.. 

Quasiparticle scattering may additionally contribute to the anisotropy in Pi{H). An 
estimate can be worked out by comparing the anisotropy in Pn{J) to that in the slope of the 
critical field /fc2(T) at T^: dHc2,c/dT ^ -7.2 T/K and dH^2,JdT ~ -4.6 T/K §. These 
values yield for the ratio of the coherence lengths 7 = C,a/^c — 0.64, which corresponds 
to an effective mass ratio of ml/ml = 7^ ^ 0.41. The anisotropy in Pn{J) amounts to a 
ratio Pn,c/Pn,a — 0.33 - 0.37. This value includes the anisotropies in effective masses and 
in quasiparticle relaxation times. Comparing the two estimates we conclude that not much 
can be allowed for the anisotropy in quasiparticle scattering. 

Finally, one more source for anisotropy is a different vortex core structure in the two 



orientations. The 2-dim Ginzburg-Landau calculations in Ref. |Ty] suggest the existence of 
several competing core structures in the axial direction H \\ c and thus the perpendicular 
direction is likely to be again different. In Fig. ^ the large Pf/pn measured with H \\ c 
could then be explained by assuming rT, = r*j_ in Eq. (H) and b?, = Sb^ = S^a^c rather than 
6| = ^l, as was done above. If we ascribe in this way the entire observed anisotropy to an 
increase in the effective core radius in the axial orientation and take 6|| =^, we estimate the 
ratio of the fourth order Ginzburg-Landau coefficients to be /52//?i = ^l/^^ ^ 0.2. This is 
in fair agreement with the independent determination (32/ (3i ~ 0.2-0.5 extracted from the 
specific-heat jump at Tc @,|]. 



Such an interpretation assumes that the effective core radius in the axial orientation 
acquires the value ^ ~ -\/3|alc ~ 2^a. This is possible if the singular "hard vortex core" has 
reduced rotational symmetry, as is the case in ^He-B at low temperatures 
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A second 

possibility is that the singular hard core lies embedded within a larger "soft-core" which 
has an effective radius comparable to ^ |Tl|. In summary, our measurements provide new 
indication for the presence of unconventional vortex structures in UPta. 

We acknowledge instructive discussions with N. Kopnin and J. Flouquet. N.L-E. and 
B.P. thank the LTL for hospitality. This work was funded by the EU Training and Mobility 
of Researchers Programme (contract no. ERBFMGETC980122). 
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FIG. 1. Spectra of the apparent penetration depth Aac(/) = A' + iA", showing as a function 
of frequency the crossover from pinning-doniinated to free-flow response (ie. the high frequency 
plateau A" ~ A' = (5f/2). The solid curves are fits to Eq. (|l|) with Ls, Ac, and pi as adjustable 
parameters. To illustrate their relative influence on the fitting, the dashed curves represent pure 
bulk-pinning (1/Ls = 0) and the dash-dotted curves pure surface-pinning (1/Ac = 0). 



FIG. 2. Field dependence of the bulk and surface pinning strengths, 1/Ac (full symbols) and 
1/Ls (open symbols) in two crystal orientations. The different symbols denote measurements 
at different temperatures below 0.33 K. The data are derived from fits to the measured Aac(/) 
spectra, as shown in Fig. |^. Bulk pinning 1/Ac vanishes at Hc2, while the surface contribution 
1/Ls decreases more rapidly and vanishes just below the B ^> C transition. The surface-pinning 
fraction is shown in gray-scale units in the inset, together with the UPta phase diagram as a 
function of field and temperature. 



FIG. 3. Flux-flow resistivity Pf{H), as deduced by fitting the measured Aac(/) (Fig. Q) to 
Eq. ([l|). It is seen that Pi{H) becomes temperature independent below 0.3 K as a function of 
H/Hc2- The results differ when measured with H \\ c (sample B3b) or iJ" _L c (sample B22). At 
low fields this anisotropy is quantified by the initial slopes ro(T) in Eq. (|2|) (dashed lines), yielding 
ro^ = 1.6ib0.15 and roii =4.7ib0.3. The inset shows unsealed data at 0.15K such that Hc2 coincides 
in the two directions. The overshoot in pf , just below Hc2, is an artifact due to the onset of flux 
penetration at the edges of the sample. In the main frame as well as in Fig. ^ such data points 
have been omitted. 



